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Abstract. It has been suggested recently that '^-exponential' distributions which form the 
basis of Tsallis' non-extensive thermostatistical formalism may be viewed as mixtures of ex- 
ponential (Gibbs) distributions characterized by a fluctuating inverse temperature. In this pa- 
per, we revisit this idea in connection with a detailed microscopic calculation of the energy 
and temperature fluctuations present in a finite vessel of perfect gas thermally coupled to a 
heat bath. We find that the probability density related to the inverse temperature of the gas 
has a form similar to a \ 2 density, and that the 'mixed' Gibbs distribution inferred from this 
density is non-Gibbsian. These findings are compared with those obtained by a number of 
researchers who worked on mixtures of Gibbsian distributions in the context of velocity differ- 
ence measurements in turbulent fluids as well as secondaries distributions in nuclear scattering 
experiments. 



1. Introduction 

Most if not all textbooks on thermodynamics and statistical physics define temperature as 
being a quantity which, contrary to other thermodynamic observables like energy or pres- 
sure, does not admit fluctuations. Because of that, it is somewhat surprising to see papers 
with the expression 'temperature fluctuations' in their titles appearing from time to time in 
serious scientific journals on subjects as various as particle physics and fluid dynamics (see, 
e.g., [[j]|2|,[3|,f|]). Indeed, how can the temperature of a system, however small, fluctuate if one 
defines it 'as equal to the temperature of a very large heat reservoir with which the system is 
in equilibrium and in thermal contact' [^]? Also, in the case of the reservoir, how can tem- 
perature be a fluctuating parameter if its definition requires one to assume the thermodynamic 
limit, i.e., to assume that the system acting as a reservoir is composed of an infinite number 
of particles or degrees of freedom? Presumably, the thermodynamic limit should rule out 
any fluctuations of thermodynamic quantities like the mean energy or the pressure, so that if 
temperature is related to these quantities, how can it fluctuate? 

The solution to this conundrum is quite simple. First, the standard definition of tempera- 
ture found in textbooks is too restrictive: there is not one but many definitions of temperature 
and of quantities analogous to temperature, as well as many physical (non-equilibrium) situ- 
ations in the context of which these different definitions admit fluctuations [^[j. Second, the 
standard definition of temperature involving the thermodynamic limit is only an idealization, 
a "purist" definition. Real physical bodies are always composed of a finite number of parti- 
cles or degrees of freedom, which means that the concept of temperature must be applicable 
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outside the idealized realm in which it is defined if experimentalists are indeed able to mea- 
sure the temperature of real bodies in real laboratories. Physically, this means also that there 
must be a threshold number of particles or degrees of freedom (more or less precisely de- 
fined) above which a system can be measured or "felt" to have a temperature [Q]. Well above 
this number, temperature is assured to be defined and likely to be constant, while close to 
this number it may be well-defined, but may change in time or vary in space; that is, it may 
fluctuate ! 

All the studies concerned with temperature fluctuations exploit one of the above 'indents' 
to the standard definition of temperature. That is, they either consider alternatives to the stan- 
dard definition of temperature which do admit fluctuations or apply the 'thermodynamic -limit' 
definition of temperature in situations where the thermodynamic limit is assumed to be "ef- 
fectively" reached without being reached "formally", so to speak. Our aim in this paper is 
to review a number of the these alternative definitions and situations, and to dispel, in doing 
so, some of the misunderstanding and misconceptions surrounding the notion of temperature 
fluctuations. We will be particularly interested in giving a detailed calculation of temperature 
fluctuations present in a system which is commonly thought to be at constant temperature, 
namely a system composed of a finite number of independent particles (basically a finite 
volume of perfect gas) thermally coupled to an infinite-size heat reservoir at constant temper- 
ature. In the following, we will show that if instead of defining the temperature of the particle 
system simply as being equal to the temperature of the reservoir we apply the statistical defini- 
tion of temperature to the finite-size system of particles (provided that the number of particles 
is sufficiently large), then we must come to the conclusion that the system's temperature is 
fluctuating, just as its internal energy is fluctuating because of the thermal coupling with the 
heat reservoir. The temperature of the particle system, in this case, is precisely related to its 
internal energy, and can be seen as a 'microcanonical temperature' associated with the micro- 
scopic configurations of a system whose internal energy is held fixed during a period of time 
shorter than the energy fluctuations time scale. 

Our motivation for studying the temperature fluctuations of a system of particles in a 
canonical ensemble setting, and for presenting moreover this study in a book about non- 
extensive statistical mechanics is threefold. First, a system composed of a bunch of inde- 
pendent particles coupled to a heat bath is one of the few thermodynamic systems for which 
the probability density describing the temperature fluctuations can be calculated directly us- 
ing 'first-principle' or 'microscopic' arguments. Second, for this specific system, the prob- 
ability density of the temperature happens to be very similar to a class of x 2 densities of 
temperature fluctuations recently introduced by Wilk and Wlodarczyk [|| S, [To||, a s well as 
by Beck [|Tl|,|l2|,[l3|], in the context of non-extensive statistical mechanics Jl4|, |15[|. Finally, 
the models of non-extensive behavior proposed by these authors are all based on the idea of 
'mixed equilibrium states,' i.e., near-equilibrium states of systems characterized by fluctuat- 
ing temperatures. In the context of the present study, this idea, as we will see, arises very 
naturally. 

2. Phenomenology of temperature fluctuations 

One can imagine many different systems exhibiting temperature fluctuations. The common 
characteristic of all of these systems is that they are non-equilibrium systems. Below, we 
list and briefly comment four systems or, more precisely, four generic situations for which 
temperature can be defined and be thought to fluctuate. The list is far from being exhaustive: 
the first three situations are presented to give an idea of the physical phenomena involving 
temperature fluctuations which have been discussed from the point of view of non-extensive 



Temperature fluctuations 



3 



statistical mechanics recently. The fourth and last case of the list, the particles and heat bath 
system, is the focus of this paper (see Section 

• Temperature fluctuations in a gas. A system with fluctuations of temperature 'spread' 
over space can be constructed simply in the following way. Take a vessel of gas, and di- 
vide it in some number of compartments thermally insulated from one another. Bring the 
content of each compartment at different temperatures, and then remove the insulating 
partition. From the moment where the partition is removed, a process of temperature re- 
laxation will take place, whereby the particles forming the gas will collide and exchange 
energy until a state of uniform temperature is achieved. 

The details of the relaxation process are quite complicated at the microscopic level, and 
depend on the nature and properties of the gas considered. But, at the macroscopic 
level, the net result of this experiment is simply described: between the time where the 
partition is removed and the time where the gas' temperature is completely uniform, the 
temperature field of the gas will vary in space as well as in time. Thus, as a whole, the 
vessel of gas can be said to be in a state of fluctuating temperature. This is admittedly an 
expletive way to say that the temperature is not homogeneous in space, but the expression 
is nonetheless correct and widely used (e.g., when referring to the spatial temperature 
fluctuations of the cosmic background radiation). 

Experimentally, there are various ways by which one can reconstruct the temperature 
field of the gas, apart from plunging a thermometer into it at different places. A simple 
method (conceptually, not experimentally) consists in measuring the momenta (along a 
fixed direction) of many particles of the gas at one point in space or, equivalently, sam- 
ple the momentum of a single particle over some period of time, and construct from the 
measurements a histogram of the number of particles L(x) having a momentum value 
between x and x + Ax. {Ax is the coarse-graining scale at which two particles are 
considered to have different momentum values.) If the sample of measurements is large 
enough, then it is expected that the form of L(x) should approximately be Gaussian, as 
predicted by Maxwell and Boltzmann, with a variance proportional to the temperature T. 
Hence, fitting L(x) with a Gibbs distribution proportional to exp(— j3x 2 ) or calculating 
its variance or its half-width all constitute operational procedures for probing the tem- 
perature T — (ksP)^ 1 of the gas. It should be noted that the accuracy of any of these 
methods for obtaining (3 depends on (i) the number of measurements used to construct 
L(x), which should be large, but not necessarily infinite!; and (ii) the assumption that 
the gas is non-interacting (perfect) or weakly interacting. These two points are necessary 
to assume that the momenta of the particles are Gaussian distributed Jl6|]. 

• Velocity temperature in turbulent fluids. It is common in turbulent flow experiments 
to define an analog of temperature by looking at the distribution L(x) of particle ve- 
locity differences in a restricted region of a fluid using anemometry or interferometry 
equipments [ |l7[ ] . Just as in the case of the gas, temperature is defined for a fluid by fit- 
ting L(x) with a Gibbs distribution of the form e~^ u ^ x \ where u{x) is the one-particle 
energy function taken to be a quadratic or a nearly quadratic function of the velocity 
variable. This defines a local inverse temperature f3 which, it is important to note, does 
not represent the physical inverse temperature of the fluid. Rather, it is a correlate of the 
local rate of energy dissipation that takes place at the microscopic level over a time scale 
known as the Kolmogorov time fl3|,|lT|], 

The Gibbsian character of L(x) and the fluctuations of the velocity temperature in space, 
related to the spatial fluctuations of the local energy dissipation rate, have been observed 



in many experiments of weakly turbulent fluids (see, e.g., 018, 19] and references cited 
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therein). However, for fluids at high Reynolds number, i.e., highly turbulent fluids, a 
totally different behavior of L(x) is observed. Indeed, recent experiments have demon- 
strated that L(x) in strong turbulence regimes is not Gibbsian; instead, it takes the form 
of a power-law which appears to be well-fitted by a so-called g-exponential function 

e-^) = [l-(l-g)^)] 1/(1 " 9) , (1) 
where (3~ 1 is a fitting parameter analogous to temperature [[nj |l8[ [l9|] . To account for 
this non-Gibbsian behavior, Beck has suggested to interpret g-exponential distributions 
as 'mixed' distributions arising from an ensemble of exponential distributions e~P u ( x ' 
parameterized by a fluctuating inverse temperature [3 flll] , |l2| , |l3| ], That is to say, if one 
assumes that what is probed in those experiments is not one velocity distribution L(x) 
characterized by a fixed temperature, but a continuum of distributions L(x) having dif- 
ferent temperatures, then what should be observed physically is an average Gibbs distri- 
bution, the average being performed over the temperature fluctuations. In this context, 
the essential point made by Beck (see [ |TTj ] for the details) is that, if the probability den- 
sity /(/3) ruling the temperature fluctuations has the following form: 



i 



1 i 

pi- 1 exp 



(q ~ 1)00 



(2) 



where (3 > and q > 1, then the mixed distribution obtained by averaging the Gibbs 
kernel e - ^"^ with /(/3) is q -exponential. Indeed, one can readily verify that 



e 



o u(x) = / e -M*lf(p)dp (3) 







using the above variant of the \ 2 or gamma density [|7|| for f(j3). This integral rep- 
resentation of the ^-exponential function is sometimes referred to as Hilhorst's for- 
mula 13 fl). 

Nuclear collision temperature. The basic idea involved in the definition of temperature 
in nuclear scattering experiments is to consider the set of particles produced during a 
collision (called the products) as forming a gas of particles which, at a first level of ap- 
proximation, can be treated as being non-interacting (perfect gas approximation). From 
this point of view, a concept of 'collision temperature' is defined essentially in the same 
way that temperature was defined for turbulent fluids except that the precise physical 
property to look at in scattering experiments is not the shape of the momenta distribution 
itself, but the so-called exponential dependence of the distribution of secondaries with 
respect to transverse momentum [[![ |l3]] . 

Since the number of particles probed during one scattering experiment is never very large 
(~ 10 — 1000), one must sometimes collect the momenta of particles over many scat- 
tering experiments before the exponential shape of the secondaries distribution reveals 
itself. However, this is not always the case: in heavy-ion experiments at very high energy, 
for example, it is often observed that a single event, i.e., only one scattering experiment 
is sufficient for a thermostatistical analysis to be effective Also, what is often seen is 
that scattering events of same nature repeated over time yield different collision temper- 
atures, making obvious that temperature is a fluctuating parameter. 
Observations of 'non-extensive' behavior in relation to this thermodynamic picture of 
scattering experiments have been reported so far on two different fronts. The first is 
related to the distribution of secondaries, and, more precisely, to observed deviations of 
this distribution from its expected exponential form. Due to the limited space available 
here, we will not discuss this case as it is quite involved. Let us only mention that Wilk 
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and Wlodarczyk have advanced in [^p a 'mixed exponential distribution' model of these 
deviations analogous to the one suggested by Beck. 

The second case of 'non-extensive' behavior concerns the absorption of cosmic ray par- 
ticles in lead chambers []8|,P,|lC[], This case was also studied by Wilk and Wlodarczyk 
who suggested for its explanation yet another variant of the \ 2 temperature fluctuations 
model (actually before Beck applied similar ideas to the study of turbulent fluids). The 
physics explained by their model is the following. The number N of hadronic particles 
absorbed in lead chambers is usually measured to be distributed as a function of the depth 
I according to 

- oc e-'A (4) 

where A is the mean free path parameter or mean penetration depth (an analog of tem- 
perature). This exponential distribution is, at least, what is observed at small penetration 
depths (~ 60 cm of lead); beyond that, what is observed is that dN/dl changes to a 
power-law which can be fitted by a g-exponential with q ~ 1.3. To account for this 
crossover, Wilk et al. simply conjectured that the A parameter characterizing the long 
flying components (i.e., the deep penetration events) is subject to fluctuations, and, thus, 
that the q-exponential penetration profiles observed experimentally for these components 
are mixtures of exponential distributions. By assuming that the probability density of A 
is a x 2 density, they were effectively able to reproduce the non-exponential distributions 
measured in laboratories [|[^J^]. 

System coupled to a heat bath. Our last example in the panorama of thermodynamic 
systems characterized by temperature fluctuations is the prototypical system defining the 
canonical ensemble: that is, a small system S in thermal contact with a larger system 
R acting as a heat reservoir. Following the standard textbook definition of the canonical 
ensemble, one should say that the temperature of system S at equilibrium is constant, 
and is equal to the temperature of system R; after all, this is how thermal equilibrium 
is defined. However, such a statement does not do justice to one important property of 
S which is that the energy density of S fluctuates (because of its finiteness) while the 
energy density of R does not (by definition of a heat bath). 

To make this statement more precise, suppose that S consists of n independent particles 
whose energy density or mean total energy is given by 

1 " 

U n = -Y j u i . (5) 

i=i 

Since the particles are coupled to R, the u^'s above are random variables, which means 
that U n is also a random variable. Moreover, observe that U n , for any finite n, has 
a non-negligible probability to assume many different values because, in this case, the 
probability density <?„(u) of U n is not a Dirac-delta function. The Dirac density, for- 
mally, is only a limiting density which "attracts" g n as n — > oo . (This basically follows 
from the law of large numbers.) Thus, if we can associate an inverse temperature fj(u) 
to all energy states such that U n = u, e.g., by applying the equipartition theorem or by 
fitting a distribution of energy levels with a Gibbs distribution as described earlier, then 
we must conclude that there are different values of (3 effectively realized 'in' or 'by' the 
particle system, so to speak. That is to say, the probability density f n ((3) for f3, obtained 
from g n (u) by a change of variables u — > f3(u), cannot be a Dirac-delta function if g n (u) 
is not itself a delta function. It is to be expected that f n {0) — > 5((3 — /?o), where {3q is the 
inverse temperature of heat bath, only in the thermodynamic limit where n — > oo. These 
points are discussed in more mathematical details in the next sections. 
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3. Energy and temperature fluctuations in the canonical ensemble 

Our analysis of energy and temperature fluctuations of a system coupled to a heat bath will 
be presented in the context of the following model. Let a vessel of gas containing n inde- 
pendent (classical) particles be thermally coupled to a heat reservoir characterized by a fixed 
inverse temperature (3q. The state of each particle is represented by a random variable Xi, 
i = 1,2, ... ,n, to which is associated a (one-particle) energy upQ). The set of outcomes 
of each of the Xi's (the one -particle state space) is denoted by X. With these notations, the 
energy density or mean energy of the gas is written as 

n 

U n {x n ) = -Y^uixi) = V L„(a;)u(a;), (6) 

i=l x<£X 

where x n = x±,X2, ■ ■ ■ , x n is the joint state of the system, i.e., the state of the system as a 
whole. Note that in the above expression we have defined L n (x) as the relative number of 
particles which are in state x, i.e., as 

T , \ # (particles : X, = x) 

L n [x) = . (7) 

n 

It should be noted that the vector L n is nothing but the histogram of one-particle states referred 
to as previously when we discussed temperature fluctuations. Indeed, in the case where x 
represents a momentum variable, the quantity nL{x) precisely counts the number of particles 
having a momentum value equal to x. (We assume throughout that the X^s are discrete 
random variables; the continuous case can be treated with minor modifications.) 

Now, owing to the fact that the gas is treated in the canonical ensemble, in the sense that 
it is coupled to a heat bath, we have 

e —/3 nU n (x n ) 

Pn(x U ) = P n {X 1 = Xl ,X 2 =X 2 ,...,X n = X n )= (8) 

AiiPo) 

as the joint probability distribution over the states x n , where 

Z n (fio)= E e-*""«W (9) 

is the n-particle partition function. Of course, since all the particles are assumed to be inde- 
pendent (perfect gas assumption), as well as all individually coupled to the same heat bath, 
we can also write 

e — 0ou(xi) e -l3 u(x 2 ) p—Pou{xn) 

P n (x n ) = p(xi)p(x 2 ) ■ • ■ p(x n ) = -—(10) 

with Z(/3q) = Zx(Po) (one-particle partition function). These equations make obvious the 
fact that what we are dealing with is a system of independent and identically distributed (IID) 
random variables. 

The first quantity that we are interested to calculate at this point is the probability distri- 
bution or probability density g n (u) associated with the outcomes U n = u. A priori, finding an 
exact expression for g n (u) is not an easy task, even though U n is the simplest sum of random 
variables that one can imagine, i.e., one involving IID random variables. Fortunately, there 
exists a general method by which one can obtain a very accurate approximation of g n (u) for 
n ^> 1 without too much efforts. This method is based on the theory of large deviations [|23|], 
and proceeds by observing that probability densities of normalized sums of IID random vari- 
ables, such as the one defining U n , satisfy two basic properties: (i) they decay exponentially 
with the number n of random variables involved; and (ii) the rate of decay is a function of the 
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value (outcome) of the sum alone. In the present context, this means specifically that g n (u) 
has the form 

g n (u)^e- nD ^. (11) 

The sign 'x' above is there to emphasize that the large deviation approximation of the density 
g n (u) is 'exponentially tight' with n, i.e., that it is exact up to 0(n~ 1 Inn) marginal correc- 
tions to the rate of decay D(u). This rate of decay or rate function is itself calculated as the 
Legendre transform of the quantity 

AO) = In E[e ku{x) ] = In ^ p(x)e ku ^ (12) 
xex 

which is the cumulant generating function of the probability distribution p(x) associated with 
the IID random variables. The result of this transform is 

D(u)=uk(u)-X(k(u)), (13) 

k(u) being the solution of 

dX(k) 



dk 



= u. (14) 

k(u) 

A proof of this result can be found in [^3[ [24[ |25[ [26|] (see also the notes contained in [ ^3| ] for 
a historical account of the developments of the theory of large deviations together with a list 
of the founding papers of this theory). 

Physicists who are not familiar with the formalism of large deviations will probably look 
at the above formulae for calculating g n (u) as being quite formal if not fancy (in a pejorative 
way). For them, we offer the following alternative derivation of g n (u). Consider the density 
Sl n {u) of states x n having the same energy U n (x n ) = u. Following the thermostatistics of 
Gibbs and Boltzmann, this density of states must be an exponential function of n taking the 
form 

n n (u)-e nH ^ u \ (15) 

where H{u) is the entropy of the system at energy density u. As is well-known, the function 
H(u) is also obtained by a Legendre transform, this time involving the logarithm of the one- 
particle partition function or free energy. Now, using the above approximation for O n (u), and 
the fact that all states x n such that U n (x n ) = u have the same probability 

e —/3 nu 

P n { x ":U n {x n )=u) = -—, (16) 

we can write 

g n (u) = Q n (u)P n (x n : U n (x n ) = u) X e -nWo+lnZ( M -H(u)] 

Thus, we arrive at 

D{u) =u0 o + bxZ{/3 o ) -H{u). (18) 

One can verify that the above expression for the rate function is totally equivalent to the one 
found in the context of large deviation theory. Both expressions are, in fact, related by the 
transformation (3{u) — /3 — k(u), where 

«.) = ^ (.9) 
du 
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is the usual thermostatistical definition of the inverse temperature. The proof of this equiva- 
lence result follows, essentially, by noting that 

A(fc) = In V 6 P ° U(X) e ku ^ = In Z(/3 - k) - In Z((3 ), (20) 

and by using the familiar expression H(u) = u(3(u) + In Z((3(u)) for the entropy. (The 
complete verification of the result is left as an exercise to the reader.) 

Let us now turn to the matter of defining an inverse temperature j3 for our system of IID 
particles, and to the complement matter of inferring the probability density f n (P)- Following 
our discussion of temperature fluctuations, it should be expected that there are many ways by 
which one can assign a temperature to the microcanonical set of states defined by 

M n (u) = {x n : U n (x n ) = u}. (21) 

Also, it is to be expected that one definition of temperature may not necessarily coincide with 
another in the case of finite-size (n < oo) systems. We illustrate this possibility by comparing 
below four different definitions or 'flavors' of temperature. 

• Derivative of entropy or free energy. An obvious way to associate a temperature to the 
states in M(u) is to take the energy derivative of the microcanonical entropy H(u) as in 
Eq.([l9|). Equivalently, one can solve the equation 

df3~ = (22) 

for /3, or compute the function k(u) from Eq.(p4|) and use the relation (3{u) = /3 — k(u). 
The inverse temperature obtained by any of these methods will be denoted by flth(u) to 
emphasize that it is based on intensive thermodynamic potentials which do not depend 
on n. 

• Derivative of the density of state. A slightly different definition of inverse tempera- 
ture is obtained by taking the 'logarithmic derivative' of f2 n (it) with respect to the total 
energy nu 

\l n (u) d{nu) dynu) 

in lieu of the derivative of the entropy exponent as in Eq.([l9| ). This defines another 
inverse temperature Pn(u) which differs from (3th{u) by a term of order (^(n™ 1 Inn) 
which vanishes as n — > oo. 

• Gibbsian distribution of states. An inverse temperature /3l(u) can be defined from a 
phenomenological point of view by fitting a given distribution of states L n (x) of mean 
energy U n = u with a Gibbs distribution of the form 

-P L {u)u(x) 

L u {x) = ■ (24) 

Z{/3l{u)) 

We have described this definition of temperature earlier (see Section ||), and have noted 
that it is accurate when n is large. To be more precise, it is accurate in a probabilistic 
sense because, in theory, there is always a possibility that non-Gibbsian distributions 
L n (x) of mean energy U„ — u can be observed. However, the probability associated 
with such a possibility is very small and vanishes rapidly as n — > oo. To see why, let us 
consider all the states x n and their corresponding distributions L n present in the energy 
'box' M{u). What we want to show is that the probability P n (L u ) that L u is observed 
in M{u) is overwhelmingly large compared to the probability P n (L) to observe any 
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other distribution L ^ L u . To show this, we use another result of the theory of large 
deviations [|24|,|2^,p^] which states that 

P n(L) _ e - n [H{V)-H{L)] _ e -nAH^ ^ 



where 



P n (L u ) 

H(L) = - ^ L ( x ) ln L ( x ) ( 26 ) 

x<£X 

is the Boltzmann-Gibbs-Shannon entropy, and AH = H(L U ) — H(L). Using the fact 
that L u is a maximum entropy distribution under the constraint U n = u, it is easy to 
see that AH > with equality if and only if L — L u , so that P n (L)/P n (L u ) -> as 
n — > 0. Moreover, the discrepancy between the two probabilities is exponentially large 
in n. Thus, for n large it can be said that any distribution L n picked at random in M (u) 
will be such that L n ~ L u . As this holds for any M(u), this implies that any measured 
distribution related to some randomly chosen state x n with b>1 ought to be a Gibbs 
distribution or be very close to a Gibbs distribution with a probability nearly equal to 1 . 

The preceding paragraphs show that there is some arbitrariness in defining the concept 
of temperature for systems composed of a finite number of particles or degrees of freedom. 
In theory, there is some indeed; however, if n is large, then defining the temperature in any 
of the ways described above should have little effect on the actual value of the temperature 
inferred. Thus, for all practical purposes, we can assume that Pth(u) — f3n(u) — (3l{u) 
for n ^> 1. In view of what was said in Section ||, it should be noted that the particular 
approximation Pth( u ) — Pl{ u ) is of deep consequences: if we look at the distributions 
L n (x) associated to the states x n G X n , then we are likely to realize that the majority of 
these distributions, i.e., those which have the most probability to be observed, form a set of 
Gibbs distributions L u parameterized by a fluctuating inverse temperature (3{u) (from now 
on we do not distinguish between the different flavors of inverse temperature). This means 
that for 1 <C n < oo all the statistical and thermodynamic properties of our system can be 
described, in an effective manner, using an ensemble of Gibbs distributions with a fluctuating 
temperature. The probability density f n {0) ruling the inverse temperature fluctuations must, 
in this case, be given by 

du(P) 



UP) = 9nHf3)) 



(27) 



where w(/3) is the inverse function of (3(u). Using this density, one then defines a 'mixed' or 
'average' Gibbs distribution of one-particle states as follows: 

L<Vf n {p)dp = J L——fJfi)dp. (28) 

This integral is a definite integral which must be evaluated over the range of definition of (3. 
Equivalently, the average can be taken over the energy coordinate: 

L(x) = J L u g n (u)du = J 6 9n{u)du. (29) 

In the above equation, be sure to distinguish the energy function u(x) from the value u of the 
mean energy U n . Also note the slight difference between these mixed distributions and those 
proposed by Wilk et al. and Beck: in our version of mixed distributions, we take the average 
over the Gibbs factor e~^ u ^ normalized by the partition function which is itself a function 
of (3 (compare Eqs.(| ) and @). 
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4. The case of the perfect gas 

As an application of the large deviation formalism, we carry out in this section the complete 
calculation of g n (u) and /„(/?) foiu(x) = x 2 /2. By using this form of energy, we assume that 
the particles composing the gas have a unit mass, and that their momentum Xi, i = 1, 2, . . . , n, 
is confined to one dimension (X is the real line extending from — oo to +oo). We also abstract 
out the position of the particles from the analysis, since the mean energy U n of the n particles 
does not depend on the position degree of freedom. 

To find g n (u), we first calculate the rate function D(u) using the Legendre transform 
method. The cumulant generating function associated with the quadratic energy function is 
calculated to be 

/°° -/3 x 2 /2 i a 

e — — e WV (30 ) 
-oo z iPo) 2 Po-k 

From this equation, we find the 'translated' inverse temperature k(u) by solving 



d\(k) 



dk 
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H u) 2 0o - k(u) 
The solution is k(u) = (3q — (2u) _1 , so that 



(31) 



D(u) = uk{u) - X(k(u)) = u/3 - - - - In 2f3 u. (32) 

Thus, 

fcf^""^"*. (33) 

This form of density is a variant of the \ 2 or gamma density mentioned previously with n 
as the number of degrees of freedom p7|]. Note that this density for the mean energy can be 
derived directly by noting that U n , for u(x) = x 2 /2, is a normalized sum of squares of n 
Gaussian random variables. In statistics, this is usually how the x 2 density is introduced [E7I], 
At this point, the density f n (f3) describing the fluctuations of (3 is readily deduced from 
the expression of g n (u) found above by calculating the physical inverse temperature (3{u). To 
this end, we can use the fact that j3(u) = /?o — k(u) or use the equipartition theorem to find 
in both cases that f3(u) = (2m)" 1 . Hence, following Eq.@, f n (/3) must have the form 

Normalizing this expression for (3 € [0, oo) yields 



= r% {^J 0- n/2 ~^~^ (35) 

as the large deviation (Id) approximation of f n (P)- A plot of this density for two values of n 
(10 and 100) is shown in Fig. [j] with (3q = 1. The plot corresponding to n = 10 has no real 
physical significance, since the large deviation approximation is not expected to be effective 
in this case. However, it is presented to illustrate the skewness (to the right) of f n (j3) which 
disappears as n — » oo. The maximum value of / n (/3) is given by /3 max = Pan/ (4 + n). As 
expected, /„(/3) converges (in a uniform sense) to the thermodynamic-limit density foo{P) = 
8{(3 — Po) when n — * oo; this is partially seen by looking at Fig. [|. By virtue of the law of 
large numbers, g n (u) must also converge in the same limit to a 6 density taking this time the 
form goo(u) = 5(u — w(A))) where 

poo -l3 x 2 /2 2 i 

<M = E H X) ] =J_ ao ^ WYdx = -. (36) 
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Figure 1. Probability densities /^ d (/3) characterizing the /3 fluctuations of n = 10 and n = 
100 free particles thermally coupled to a heat bath with /3q = 1. Each density is defined for 
/3 > 0, and shows a maximum at /3nTl/ (4 + n). 




Figure 2. Comparison of the x 2 /3-density proposed by Beck (B) Jl 1[ and the one proposed 
in this work (T) for n = 10 (left) and n = 100 (right) particles (see text). For n = 100, the 
two densities are quasi-indistinguishable. The maximum value of the j3 -density, in the case of 
Beck, is located at /3o(n — 2)/n. 



We now come to the main point of our study which is to compare the f n (P) density ob- 
tained here and the x 2 probability density of Eq.(g) which was 'postulated' by Wilk etal. and 
Beck in their studies of mixed distributions (see Section ||). To establish this comparison, we 
present in Fig. || two plots of fjf(P) for two different values of n and a variant of the \ 2 P 
-density proposed by Beck 

1 / \ n /2 
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which results from identifying l/(q — 1) in Eq.(^J) with n/2 [|ll[]. The plots are presented 
again for n = 10 and n = 100. A rapid inspection of the expressions of and / r f reveals 
that these densities are not at all the same. The density can in fact be viewed as emerging 
from the sum of the squares of n Gaussian random variables, whereas, in our case, the sum 
of squared Gaussian random variables arises as the mean energy, and so as modulo some 
constant. This explains why performing the change of variables [3 — > /3 _1 in f~f yields /,f 
modulo some constant and a Jacobian term arising from the change of variables. In spite of 
this important difference, the second plot of Fig. ^ shows that both densities are remarkably 
similar as n gets large. This, at first, does not seem surprising as both densities converge to the 
delta density foc(P) in the thermodynamic limit n — > oo. However, it is to be noted that each 
of them gives totally different mixed distributions when they are used in Eq.(^8|). Indeed, it 
can be shown [ ^8| ] that the mixed distribution associated with f,f has the asymptotic form 

L l *(x) ~ e-W (38) 

for \x\ ^> 1, instead of 

L B (x)(xe~ x2 , (39) 

where q = 1 — 2/n. Both of these results should be compared with the (pure) Gibbsian 
distribution 

L G {x) cx e~ x2 (40) 

which is the limiting distribution of L~* and L B in the thermodynamic limit (n — > oo). 

The above scaling relationships clearly indicate that choosing between /„ or /jf 
has a dramatic consequence on the functional form of the mixed distribution calculated even 
for n ^> 1 . Does that imply that our model of temperature fluctuations cannot serve as a model 
of the non-Gibbsian distributions which have been observed in turbulent fluid experiments as 
well as in nuclear scattering experiments? The answer is not as straightforward as one would 
think. First, it is not at all clear that turbulent fluids can actually be treated in the canonical 
ensemble and/or that the perfect gas assumption is a valid approximation in this case. These 
points call for further justifications. Second, the extreme events |a;| ^S> 1 needed to validate 
either one of two temperature densities compared in this work are often very difficult to detect 
experimentally in a reliable way. Surely, additional calculations and experimental data would 
be welcome in order to test the validity of our approach, and to confront it with that of the 
authors mentioned in the present study. This seems to be especially true for nuclear scattering 
experiments which are usually thought to fit perfectly well into the canonical ensemble picture 



5. Concluding remarks 

Our treatment of the energy and temperature fluctuations of a system coupled to a heat bath 
has focused mainly on the perfect gas. However, it is worth noting that the large deviation 
approach presented in this paper for calculating the energy and temperature probability den- 
sities in the canonical ensemble is very general. It can be applied independently of the form 
of the energy function u(x) which defines the mean energy U n , and can also be generalized 
without too much difficulties to cases involving other forms of probability distribution for 
Pn(x n ) (e.g., q -exponential distributions). In this context, an obvious extension of our work 
could be to consider different forms for u(x), and to look at the mixed distributions which re- 
sult from the corresponding temperature fluctuations. This line of thought has been followed 



recently by Beck and Cohen [30] who derived a number of 'superstatistical' mixed distribu- 



tions (sometimes unphysical ones) by assuming different forms of temperature fluctuations. 
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Another problem could be to solve the following 'inverse problem': for which u(x) is f n (j3) 
the same \ 2 density as the one suggested by Beck? Finally, note that a large deviation cal- 
culation of g n (u) and /„(/?) can also be carried out for systems involving dependent random 
variables. Unfortunately, the calculations leading to the specific forms of g n (u) and f n {0) in 
this case are likely to be tedious. Also, the concept of mixed distribution does not generalize 
easily to the case of interacting particles because the Eqs.(^5|) and ( p6[ ) which were used to 
prove that Gibbs distributions are the only distributions likely to be observed in large systems 
are valid for sequences of IID random variables only. It is, in fact, a long-standing open prob- 
lem of large deviation theory to generalize these equations to sequences of dependent random 
variables. Solving this problem would have direct consequences in statistical physics, for it 
implies ipso facto a generalization of the maximum entropy principle to systems of interacting 
particles. 
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